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Abstract 

We show that at any temperature, the low-energy (with respect to the chem- 
ical potential) collective excitations of the transverse components of the energy- 
momentum tensor and the global U(l) current in the field theory dual to the 
planar RN-AdS4 black hole are simply those of hydrodynamics. That is, hydro- 
dynamics is applicable even at energy scales much greater than the temperature. 
It is applicable even at zero temperature. Specifically, we find that there is always 
a diffusion mode with diffusion constant proportional to the ratio of entropy den- 
sity to energy density. At low temperatures, the leading order momentum and 
temperature dependences of the dispersion relation of this mode are controlled 
by the dimension of an operator in the thermal CFTi dual to the near-horizon 
Schwarzschild-AdS2 geometry. 
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1 Introduction 

Over the last decade, it has become clear that gauge/gravity duality [TH3] is a very useful tool 
for studying the hydrodynamics of strongly-coupled field theories (see [1H6] for some reviews 
of this field). The theory of hydrodynamics [7j is formulated as a derivative expansion in 
which the macroscopic quantities describing the state of a system close to equilibrium are 
assumed to be slowly- varying with respect to a length scale / mfp , often taken to be the mean 
free path between thermal collisions. Typically, increasing the temperature T will reduce 
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/ m f p and thus increase the accuracy of hydrodynamics. A recent introduction to this topic 
can be found in [8]. 

The theory of hydrodynamics predicts that for a conformal theory in D spacetime dimen- 
sions deformed to non-zero T and with a chemical potential fi for a conserved U(l) charge, 
the two-point functions of the energy-momentum tensor T^ u and of the U(l) current J 11 ex- 
hibit three kinds of gapless excitation: a transverse diffusion mode with dispersion relation 
(where u is the frequency of the excitation and q its momentum) 

^ = -*7^ 2 + o (<7 4 ) , (1) 
a longitudinal charge diffusion mode, and a longitudinal sound mode with dispersion relation 



where e is the energy density and P the pressure of the field theory (as the underlying theory 
is conformal, these are related via e = (D — 1) P). The shear viscosity rj of the theory, which 
controls the attenuation of the sound and transverse diffusion modes, can be independently 
determined via the Kubo formula 

rt = — lim -ImGLps (u, q — 0) , (3) 

uj^O CO 

where Gq^q. denotes the retarded Greens function of a pair of operators. 

Starting with the existence of these modes has been established in a wide variety 

of holographic theories in the limit w,g<T, which is analogous to the limit described above 
where the thermal collision length is much shorter than the wavelength of the perturbation. 
Remarkably, it has been shown that for a very large class of holographic theories the viscosity 
is given by [I3HIS] 

' = (4) 

where s is the entropy density of the field theory. This result is significant as the ratio 
rj/s is similar to that measured for the quark-gluon plasma, and significantly smaller from 
that predicted by, for example, perturbative QCD. One obvious lesson from this work is 
that the transport properties of strongly-coupled theories can be significantly different from 
those of theories based on the existence of long-lived quasiparticles. The applicability of 
hydrodynamics to holographic theories can be checked to higher orders in the derivative 
expansion and in the amplitude of fluctuations using the fluid/gravity correspondence [T71 - 

Q2|. 

A question which has had comparatively little study in the gauge/gravity duality litera- 
ture is what happens when the mean free path between thermal collisions becomes very long 
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- when u,q ^> T, or in the extreme limit when T = 0. One could argue that hydrodynamics 

- a derivative expansion valid on length scales much longer than Z m f p - should be valid even 
at T = provided that we consider perturbations with uj, q <C ji- In other words, that at 
T = hydrodynamics will continue to be valid since a large /i will result in a small mean 
free path Z m f p ~ and thus a well-controlled derivative expansion. This viewpoint was 
emphasised in [20] : assuming that rj/s = 1/Att at T = 0, one can define Z m f p to be equal to 
the diffusion constant in ([T]) and show that it is finite at T = for a certain class of large, 
charged black holes in global AdS spaces. 

However, this is at odds with our intuition that turning on a large T or a large \i 
are qualitatively different deformations of a field theory. Moreover, in the best-understood 
quasparticle-based theory of this kind - Fermi liquid theory - the hydrodynamics outlined 
above is not valid at sufficiently low temperatures^ For example, the speed of sound at 
zero temperature is not equal to its hydrodynamic value \JdP/de and the decay rate of 
sound — Im (u) oc q 2 + T 2 has an important, g-independent contribution [2"TH2~4"] . Further- 
more, it has been shown in various holographic probe D-brane theories that hydrodynamics 
is not an accurate description at sufficiently low temperatures, as the longitudinal charge 
diffusion mode expected from the hydrodynamics of a conserved U(l) current is replaced by 
propagating sound modes [25T436] . 

There have been a few studies of the low-energy bosonic excitations outside of the usual 
hydrodynamic limit u,q < T for non-probe holographic theories with a non-zero fj,, including 
[37H4Tj . The main technical difference of these calculations, compared to those in the usual 
hydrodynamic limit, is that as T — > the location of the black hole horizon becomes an 
irregular singular point of the equations of motion for the excitations of the bulk fields. 

Additionally, it has been shown that if one defines the 'viscosity' at T = formally using 
the Kubo formula <^), then for a small subset of holographic theories (which all have s ^ at 
T = 0) the relation rj = s/4tt continues to hold even at T = [I2]-FJ6]. This subset includes 
the planar AdS-Reissner-Nordstrom black hole. We stress here that these examples are not 
covered by the general proof of [16] because of the different horizon structure in the T = 
limit mentioned above. Physically, the two situations are completely distinct as the u — > 
limit in the definition of rj corresponds to w < T, /i at any non-zero T, and T <^ u> ^ fi 
when T = 0. Henceforth we will abuse notation and denote both the usual viscosity, and 
its zero temperature analogue, by the name 'viscosity' and the symbol rj since they are both 
equal to s/Att for the field theory that we investigate. 

In this paper, we study the retarded Greens functions of the transverse components of 
T^ u and J M in the field theory dual to the planar AdS-Reissner-Nordstrom black hole in 
(3+l)-dimensions (RN-AdSzi). We analytically compute these Greens functions at leading 
order in small u, q under the assumption that these quantities are small compared to either 



1 It breaks down when uj > l m ^ p ~ T 2 /fi. 
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T or fi (or both). When u,q <C T, we recover the expected hydrodynamic diffusion mode 
with a dispersion relation given by ([Q) and (]4]), in agreement with the results of [UJHB]- More 
interestingly, we show that this diffusion mode exists with a dispersion relation given by ([I]) 
and 01]) even when u > T. In fact, even in the extremal limit T = 0, we find that there exists 
a diffusive mode described by the hydrodynamic results ([I]) and (jlj). To our knowledge, this 
is the first analytic calculation of the bosonic excitations outside of the usual hydrodynamic 
limit w < T in a holographic theory at non-zero density with a dynamical metric. 

Outside of the usual hydrodynamic range u <C T, our analytic calculation of the dis- 
persion relation of the diffusion mode relies on the existence of a (Schwarzschild-)AdS2 
near-horizon geometry. In the spirit of the semi-holographic work of [To ^B^HoT] , we find that 
whereas the existence of low-energy modes (in our case diffusion) is a property of the full 
spacetime geometry, the decay rate of these modes is controlled by the near-horizon geome- 
try. In our case, the leading power of the momentum in the dispersion relation u ~ — iq a is 
controlled by the dimension of the operator dual to the gauge-invariant combination of the 
metric fluctuations in the CFTi dual to the near-horizon AdS2 geometry. Formally, we find 
that the dispersion relation of the diffusion mode at T = can be written 



-q\ (5) 



where (?J R (u) is the retarded Greens function of a scalar operator of conformal dimension 1 
in the CFTi dual to the near-horizon AdS2 region of the geometry. As this CFTi correlator 
is proportional to iu, we obtain u ~ — iq 2 . Note that unlike in the fermionic case of [49] . the 
dimension of this operator is fixed. At non-zero temperatures T < uj <C fi, (?J R (u) is replaced 
by QYt the thermal Greens function of the CFTi operator, in the above equation. This 
ensures that the dispersion relation u (q) retains the same u ~ — iq 2 form (but with a In- 
dependent diffusion constant). In particular, it does not receive any corrections which are 
dependent upon T but independent of q (at lowest order). 

The T = Greens functions of the transverse operators of this theory were analytically 
studied previously in [23112], and our results are complementary to theirs. We use different 
gauge-invariant fields than in their approach, which allow us to determine the leading-order 
momentum dependence of these fields far from the horizon. As just indicated, it is the 
knowledge of this far-from-horizon solution which allows us to demonstrate the existence of 
a diffusion mode. On the other hand, the fields used in [37] |4"2] allow one to determine the 
dependence of these fields near the horizon to a higher order in u and q than in our approach. 
This means that the scaling of the spectral functions in the limit u — > is captured more 
accurately by the results of J37JI12] (although at leading order, our results agree). 

The remainder of this paper is structured as follows. In section [2] we review the ther- 
modynamic properties of the planar RN-AdS4 solution, and manipulate the equations of 
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motion (and on-shell action) of the fluctuations of the transverse fields into a form in which 
we can solve them in the appropriate limits. In section [3] we calculate the Greens functions 
of the dual operators of these transverse fields in the limit w < T. This is an instructive 
warmup for sections H] and in which the same general procedure is used to calculate the 
Greens functions in the T = (with u <C fi) and T < u <C fi limits respectively. In section 
[6] we present results of the numerical computation of the poles of the Greens functions, and 
the spectral functions, at both zero and non-zero temperatures. These confirm the accuracy 
of our analytic results. Finally in section [7] we conclude with a discussion of some possible 
future research directions. 

2 The planar RN-AdS4 solution and its linearised ex- 
citations 

2.1 The equilibrium solution 

The planar AdS4 Reissner-Nordstrom spacetime that we study is a solution to the Einstein- 
Maxwell theory with a cosmological constant. The action of this theory is 

S=—^ [ d^x^g ( Tl + ^ - L 2 F^ u F^i + \ f d z x^/\h\ K, + counterterms, (6) 

2^4 J M \ L J K 4 JdM 

and the planar RN-AdS4 solution can be written 

ds 2 = -—^p-dt 2 + r — (dx 2 + dy 2 ) + dr 2 , 
L 1 1/ r z j(r) 

3 4 

/(r) = i-(i + g 2 )^ + g 2 ^, 

where the boundary of the spacetime is at r — > oo, the planar outer horizon is at r = tq and 
the U(l) field strength is given by F = dA. The term involving the extrinsic curvature K, is 
the Gibbons-Hawking term, and we have not written explicitly the counterterms required to 
remove the divergences from the on-shell action. These counterterms, which can be found 
for example in [37] . contribute only contact terms to the dual Greens functions and so will 
not be required for our purposes. 

This solution can be embedded into 11-dimensional supergravity, where it is realised as 
the near-horizon geometry of a stack of spinning M2-branes. It is dual to the state of the 
low-energy (2+l)-dimensional field theory on a stack of M2-branes with a non-zero density 
of diagonal U(l) R-charge [521 153] . The solution has one dimensionless parameter Q, which 
determines the ratio of the temperature of the field theory state T to the chemical potential 



A t (r) 



Qrp 
L 2 



ro 
r 



(7) 
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(with respect to the diagonal U(l) R-charge) \i 



(3-Q 2 )r Qr T 3 - Q 2 

T = T~F> > V ~ 



4ttL 2 ' L 2 ft AnQ 

This can be inverted to give Q in terms of T/fi 



Q= ./3 + 4vr 2^!_27r-. (9) 

Q can take values between (corresponding to /i = 0) and a/3 (corresponding to T = 0). 

The thermodynamic properties of the dual field theory state are controlled by the dimen- 
sionless parameter Q (or equivalently by T/fi). The energy density e, pressure P, entropy 
density s and charge density p are given by [52] 



K 2r4 ' 9' t-2 r 2 ' ^ ~ K 2 j 2 ' ^ iU ^ 

When T = 0, this state has the unusual property that its entropy density is non-zero. 
2.2 Linearised fluctuations around equilibrium 

The two-point functions of field theory operators can be computed by studying the linearised 
fluctuations of their dual fields in the gravitational theory. Without loss of generality, we 
choose the momentum of these fluctuations to flow along the x-direction and excite the 
metric - which is dual to the energy-momentum tensor T^ u of the field theory - and the 
gauge field - which is dual to the diagonal U(l) R-current J M - around their values on the 
planar RN-AdS 4 solution as follows 

9,u (r) 9lw (r) + f ^e-^+^V (r, u, q) , 

J {2n) (11) 
A,(r)^A,(r) + J ^|e"^X (r, g) . 

At linear order, the fluctuations of fields transverse to the momentum flow - h yt , h xy , h ry 
and a y - decouple from the rest, due to the y — > —y symmetry of the theory and of the 
planar RN-AdS4 solution. In this paper we are interested only in the fluctuations of these 
'transverse fields'. By varying the action (jBJ) with respect to h yt , h xy , h ry and a y respectively, 
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we obtain the following set of coupled equations of motion 



±[T*fhZ'-i q L%]+^( q hi + U hi)=0, 



(qh y t +cuh x y )=0, 



V f) h y+ r 2f ' 



0. 



0. 



(12a) 
(12b) 
(12c) 
(12d) 



where we raise and lower indices of using the background metric <^ and where a prime 
denotes a derivative with respect to r. The metric fluctuations with mixed indices are 
more convenient for numerical calculations, as their leading term near the boundary of the 
spacetime is a constant. 

Not all of these equations of motion are linearly independent - equations f)12ap and 
(I12c| together imply equation (112bfl and similarly equations (112bfl and (j 12 eft together imply 
equation f)12ap . This is a manifestation of the underlying gauge symmetries of the linear 
fluctuations which fix the relationship between h v t and h x y as we will shortly demonstrate. 
For our purposes, it is convenient to choose ur 2 A' t x (112fcj) + qr 4 f x (112<ij) + qr 2 fA' t x f|12ap . 
along with equations f!12bl) and H12cj) . as the three linearly- independent equations of motion. 
This first combination can be written 



A. 

dr 



r"fA' t ( hf + -K% ) + r 6 //'< 



L A r 2 A' t 
L A r 2 f 



+ 
+ 



f 



(u 2 - q 2 f) a y 



(13) 



where we have used the fact that r A' t (r) is independent of r and that 



A 



t j 



(14) 



for our background. This is a convenient choice because the terms outside the derivative 
bracket (which are gauge-invariant) are suppressed at sufficiently low u and q with respect to 
the terms inside the derivative bracket (which are also gauge-invariant). Thus at sufficiently 
low lo and q, this equation (along with equation (I12bj) ) can be trivially integrated to produce 
gauge-invariant first-order differential equations for the fields. 

As alluded to above, the linearised fluctuations possess a diffeomorphism gauge symmetry, 
under which the fields transform as [54] 



(15) 
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where the covariant derivatives are taken with respect to the planar RN-ACIS4 metric (J7J). The 
gauge field fluctuations also transform as a M — >■ a M — <9 M A under the bulk U (1) gauge symmetry, 
which means that the transverse gauge field a y is invariant under this transformation, as are 
all of the metric fluctuations h^ u . It will be convenient to write our theory in terms of 
gauge-invariant combinations of the fluctuations of the fundamental fields h^ v and a M . One 
can construct various different gauge-invariant combinations from the fundamental fields and 
their derivatives, and the choice we will make is to consider combinations of fields which do 
not have an r index. It is these fields whose boundary values have a clear interpretation 
as the sources of the dual field theory operators. With this restriction, the gauge-invariant 
fields are 



¥1 (r, u, q) = h\ (r, u, q) + -hi (r, u, q) 



q 



(p 2 (r,U),q) = — ay(r,u,q), (16) 



or any linear combination of them. When written in terms of these variables, the gravitational 
theory explicitly encodes the Ward identities of the operators in the dual field theory [55] 



{J r £-x,y r pty — \J ' J^tyJ^xy 



{-J' r pxy J^xy 



UJ 



2 *- r T t 2'T t !' 1 



/-iR 



where these expressions should be understood to hold up to contact terms. This will be seen 
clearly in the following sections. 

To obtain the equations of motion in terms of these variables, one simply solves equation 
(Il2cj) algebraically for h r y and then substitutes this solution into the dynamical equations 
( 1X3]) and (I12bj) . This yields the following coupled equations of motion 



d_ 

dr 
d_ 

dr 



r f ( A<Pi + -JY^2 



7-4 2 



./ 



A' 1 r °f 



0. 



r 2 f 



q 2 f 



(rVi + 4r L 2 ^ 2 ) 



L 4 

+ y^i =0, 



(18a) 
'18b) 



which will be the subject of the remainder of this paper. 

To compute the Greens functions of the dual field theory, we also need to know the on- 
shell gravitational action. This can easily be computed in terms of and a^, and after 
substituting in the solution for h r y from equation fl!2cj> . it becomes 



S 



dudq 



2k\ Jr-toc (27r) 2 



r 4 fq 2 



2L 4 (co 2 - q 2 f) 
2r 2 r 2 f 



(fx (r,-u,-q) (p[ (r,u,q) 



L 4 



ip 2 (r, — u, —q) (p' 2 (r, uj, q) + non- derivative terms 



(19) 



where a prime denotes a derivative with respect to r. The non-derivative terms in this 
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action produce contact terms in the Greens functions. These non-derivative terms cannot be 
written purely in terms of the gauge-invariant variables above, which is simply a reflection 
of the fact that the Ward identities (I17p hold up to contact terms. In the gravitational 
theory, it is not a signal of broken diffeomorphism invariance but rather that the linearised 
transformations ([151) should be modified at quadratic order. 



3 Greens functions in the usual hydro dynamic limit 

As a warmup for our low temperature calculations, in this section we will determine the 
Greens functions in the usual hydrodynamic limit u <T. As expected, we find a diffusion 
pole with a dispersion relation given by equations (pQ) and @. Our results are consistent 
with those of [H], although our method is different. Our method can easily be generalised to 
compute the Greens functions outside of this limit, as we will show in the following sections. 

To find the Greens functions, we divide the bulk spacetime into an inner region (a 
suitably-defined region near the horizon) and an outer region (a suitably-defined region 
near the boundary). These two regions overlap over a range of r called the matching re- 
gion. Firstly, we solve the equations of motion in both the inner region and outer region 
to determine <pi and <£>2 in the respective regions up to integration constants. After sub- 
stituting in the solutions <y2° uter and (j( 9° uter ; the boundary on-shell action (from which the 
Greens functions are computed) is determined in terms of the integration constants of the 
outer solutions. These integration constants should be fixed by demanding that the fields 
are ingoing at the horizon. This is implemented by imposing ingoing boundary conditions 
on the solutions in the inner region, and then demanding that the solutions <f^ eT and v 9 ° U 2 ter 
are consistent in the matching region. Having fixed the integration constants of the outer 
solutions, we can read off the Green's functions from the on-shell action. 



3.1 The inner region 

We begin by solving the equations of motion (I18al) and (I18bl) in the inner region. After 
diagonalising the two-derivative terms, we expand these equations around the horizon r = ro 
(assuming that (p\ ~ (p^) to obtain 



1 



r — r 
<P2 + <Pi [•••]+ V2 



+ ¥? 2 [•••] + <Pi 



L A uj 2 



r — To 



+ 



+ y?2 



L(3-Q 2 ) 2 r 2 (r-r ) ; 
L(3-Q 2 ) 2 r 2 (r-r ) ; 



+ 



+ <?2 [• ■ •] 

= 0, 



0. 



(20) 
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where the ellipses denote higher-order terms in this near-horizon expansion. In this region, 
the fields decouple and their equations of motion can easily be solved to give the usual result 
for fields near a horizon of non-zero temperature 



, „inner 
¥>1,2 



1,2 



r 



T & 

r {3-Q' J ) 



'1,2 



r 



aj 2 exp 



4^T 



log - - 1 

r 



+ a 12 exp 



iuj 




loe 


(Hi 


4vrT 6 





(21) 



where integration constants. Ingoing boundary conditions correspond to the choice 

<2 = 0. 

Note that in the zero temperature limit (Q — > Vo), the expansion (120|) breaks down, as 
terms which naively appear to be of a higher order contain factors of the form (3 — Q 2 ) \ 
This is simply a consequence of the horizon structure changing in the T — > limit as f(ro) 
becomes a double zero. 



3.2 The outer region 

We will now solve the equations of motion ( 1 18 aft and (|18bl) in the outer region. We define 
this region by 

u 2 L 4 q 2 L* , x 

~*p « 1. 72J « 1- ( 22 ) 

After taking these limits, the non- derivative terms in the equations of motion drop out and 
we can trivially integrate them to give 

r J \ A t<fi + ~JY^2 J = L 2 » (23a) 

r ^2, outer' , r 2 outer\ C 2 r 



— q 2 f 



(r^r er ' + ^ L 2 A^r r )=^, (23b) 



where C\ and C2 are integration constants and we have used the fact that r 2 A' t is a constant 
for the planar RN-AdS 4 background. It is then simple to decouple ^° uter ? which obeys the 
following linear, first-order differential equation 

outer' i^Louter ^ (Cl ~ C 2 ) + C 2 q 2 f 
^2 '—^JT^ = • (24) 

This equation can be formally solved by means of an integrating factor to give 



outer / x _ „„„ / / A t 



^° utcr (r) = exp / dr 



f.2 jl 



b\ + I dr exp I — / df 



f AL*A> 2 \ u 2 (d - c 2 ) + c 2 q 2 f 



f J r& f f 

(25) 
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where &i is an integration constant. Substituting this solution for <^2 Uter into equation (123bj) 
we can easily integrate to obtain <^° uter 



^outer (r) = b2+ d f 



r c 2 {u 2 -q 2 f) 4r L 2 A' , ul( . 



L 2 f Q fA' t 

where b 2 is an integration constant. Doing the integral 



^2 \ r , 



(26) 



r ,AL 4 A' 2 , 



we find that 



(r) = <^ 0) + / df 



r c 2 (u 2 - q 2 f) 



3 (1 + Q 2 ) - 4Q 2 ^ 



(27) 



L 2 f 6 fA' t 
1 - 



4r L 2 ^ ^ 4Q 2 r 



3(i + g 2 ) f 



4Q 2 r 
3(1 + Q 2 ) f 
df lu 2 (ci - c 2 ) + c 2 q 2 f 



(2f 



4Q 2 r 
3(1+Q 2 ) f 



and 



¥>2 (r) = 1 



4Q 2 r 
3(1 + Q 2 ) r 



(o) , 

f2 + 



df 



u 2 (ci - pa) + c 2 g 2 / 



1 - 



4Q 2 r 
3(1+Q 2 ) r 



(29) 



where <p± \ are the values of (fi j2 at the boundary of the planar RN-AdS4 spacetime. 
Doing the integrals for the terms proportional to q 2 , these become 



^ Utcr (r) = ri 0) + ^ ( 1 



Q 2 r 



3r 3 V (1 + Q 2 )r 



Qr 



L 2 f 6 fA' t 



4r L 2 A' t 4Q 2 r 



df 



v (Ci - c 2 ) 



(30) 



3(i + Q 2 )fyy f«fp(\ AQ2 r A 

r J J \ l 3(1+Q 2 ) f ) 



and 



(r) = 1 



4Q 2 r 
3(1 + Q 2 ) r 



(o) 

¥>2 



c 2 g 



3r 3r(l + g 2 ) (l-5^y=? 



(31) 



r6 f f M — -JSLai 

\ 3(1+Q 2 ) r y J 

We could not do the integrals exactly for the terms proportional to u 2 (except in the T = 
or fi = limits), but this will not be important for determining the Greens function at lowest 
order in u and q, as will become clear shortly. 
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3.3 Matching 

To completely determine the outer solutions and hence the Greens functions, we must fix the 
integration constants c% and c 2 such that the fields are ingoing at the horizon. As previously 
mentioned, this can be done by finding a matching region - a range of r over which both the 
inner solutions and outer solutions are valid - and comparing the solutions in that region. 
As the inner region solutions are valid near the horizon, we expand the outer region solutions 
near r = tq. For the terms proportional to u 2 , it is sufficient to expand the integrands in 
this limit and keep the leading terms to give 



Qro (3 - Q 2 ) \r 



(o) . ^Qv 



(o) 

2 



+ 



C 2 <f 



3(1 + Q 2 ) 3Qr4(l + Q2) 



+ O (u 2 c 1 ,u 2 c 2 ) 



+ O (r - r ) 



(32) 



and 



outer 
<P2 



AQ 2 \ 3 (1 + Q 2 ) u 2 ( Cl - c 2 ) /r 

Hi + Q 2 )) (s-g 2 ) 3 r 4 og U 



+ 



V f I 1 



4Q 2 



c 2 q 



+ O (u 2 c 1 ,u 2 c 2 ) 



(33) 



+ O (r — r ) . 



where we have neglected terms of order ciw 2 and c^io 2 in the coefficients of the (r — r )° 
terms. This is because we are searching for a diffusive mode with u ~ q 2 and thus the bJ 2 C\ 
and u 2 C2 terms are subleading in an expansion at low q. When r — ro tq, as in these 
solutions, our constraint on the validity of the outer region solutions ( l22p becomes 



ujL 2 r — Tn 

< < 1. 



(34) 



This requires that w<T, which is the hydrodynamic limit commonly used in holographic 
computations of Greens functions. 

The form of the inner solutions (1211) in the matching region are found by expanding them 
in the limit 

(35) 



to give 



on) ( f 

(r) = K2 + % 2 ) + « 2 - a r , 2 ) — log l- 



1 + 



(36) 



Demanding ingoing boundary conditions at the horizon (a^~ 2 = 0) fixes the relative coefficient 
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of the constant and logarithmic terms in this expansion 



VST (r) = a7 2 



Sr log U _1|H 



(37) 



To determine ci and C2 for a solution which is ingoing near the horizon, we equate the 
coefficients of the constant and logarithmic terms in the outer solutions fl32l) . fl33|) with those 
in the inner solutions (13T1) . The resulting two equations are 



iuc 2 (o) . 4Q(^2 0) 



gr 3 L 2 
iu (ci - c 2 



+ 



c 2 q 



3(l + g 2 ) 3gr4(l + g2; 

4g 2 



+ O (ciw 2 , c 2 c<; 2 ) , 



^ 0) 1 1 



c 2 q 



L 2 r 3 (3-g 2 ) r2 V 3(i + g 2 ); 3^(l + g 2 ^ 

Neglecting terms of O (u 2 ci,u 2 c 2 ), one can solve these to give 



+ O (ciw 2 , c 2 co 2 ) . 



Ci 



7-23 



L 2 g 2 



3r (l+Q 2 ) 



(0) / . L 2 (3 - g 2 ) 2 
3r (l + g 2 ) y 



g^i iu; - 



+ 



(o) 

¥>2 



3(i + g 2 ) 



Lj (g 4 -2g 2 + 9)-^(s-g 2 )) 



c 2 



1+Q 2 ) 



<^2 



IU0 



L 2 q 2 



3r (l+Q 2 ) 



(38) 



(39a) 



(39b) 



These constants, along with the solutions ( 13TJ|) and (1311) . determine the behaviour of tp\ and 
(f 2 near the boundary of the spacetime, assuming that the fields are ingoing at the horizon 
and that we are in the usual hydrodynamic limit w <T. 



3.4 Green's functions 

To determine the retarded Green's functions, we require the on-shell action for fields that 
are ingoing near the horizon. Substituting the outer solutions fj30|) and fl3~T|) into the on-shell 
action (fT9]) . we find (where the argument (to) is shorthand for (uj,q)) 



S 



dudq 



1 



2kI (27r) 2 2L 2 I loJ 



-r 2 q 2 



L 2 q 2 



3r (l+Q 2 ) 



(-a;) <pY> (to) + 



(0) 



i6g 2 



9(i + g 2 



(°) i \ (°) i \ 
^2 (-w) <^ 2 M 



+ 



4g 



3(i + g 2 ) 



(40) 
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where we have imposed the values of c\ and c 2 for ingoing modes (j39aj) and (139bl) and 
neglected contact terms. Using the relations 



^ ( W , q) = hf 0) ( w , g) + -/if) («, g) , <4° } («, g) = -af (u, q) , (41) 

g r y 

where the superscript (0) denotes the boundary value of a field, and the usual prescription 
for calculating Greens functions from holography [56j[57], we find that the retarded Greens 
functions are (up to contact terms) 

r n rig- _ 8L 2 q 2 Q 2 



(42) 



r R _ r R 2r q 2 Q 
{ ^ r T t yjy — " jvT l y 



3^(1 + ^) (zo;-^ 



3n,(l+Q 2 ) 

with the Greens functions of the other operators obeying the Ward identities ( 1T7|) due to 
our gauge-invariant bulk formulation. These results are valid only at lowest order in u and 
q (assuming that u ~ g 2 ), and are valid only up to contact terms. They agree with the 
results of up to such correct ions They accurately capture the behaviour of the spectral 
function in the neighbourhood of the diffusion peak, as we will show numerically in section 

EJ 

The Greens functions share a common diffusive pole with dispersion relation 

L 2 s 
3r (1 + Q 2 ) 4:71 (e + P) 

where the ellipses denote higher order terms in q, as one expects from the hydrodynamic 
equations flTJ and (J4J). This value of the viscosity can be checked by applying the Kubo 
formula (ED to our results for the Greens functions. 



4 Zero temperature Greens functions 

As we emphasised, the above calculation of the Greens functions is not valid when T = 0. 
This is because the horizon structure is qualitatively different in this limit - f(r) has a double 
zero when T = 0, rather than the single zero it has when T ^ 0. However, at T = one 
can use a similar procedure to that just outlined provided that the inner region is suitably 
defined. For the extremal background there is a natural candidate for such an inner region 
- the AdS 2 x M 2 near-horizon geometry. 

2 Although the correlators involving J v appear superficially different from those in [?5], these differences 
are just due to contact terms. Both our calculation and the calculation in [?5] neglect certain contact terms, 
but not the same set of contact terms (since they are done by different methods). 
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4.1 The inner AdS2 region 

Following similar calculations in [371 [381 HD 1121 EH] ; the near- horizon A0IS2 x R 2 geometry is 
most easily seen by changing co-ordinates to 



c 



r — ro 



(44) 



and then expanding linear combinations of the equations of motion (11 8 aft and (118bj) as power 
series at small u and q (assuming that <px ~ <£> 2 ) to give 



(0 + <A (0 



+ ^ (C) [...] + n (C) 



(0 + <p[ (0 [•••] + ^ 2 (0 



+ ¥>2 (C) 



36C 4 

1 

36C 1 



2 



+ ^ (0 [• • •] 
0. 



0. 



(45) 



where the ellipses denote terms with positive powers of u, q. In these limits, the equations 
of motion decouple and after redefining the radial variable £ — >■ £/ (w£ 2 ) to remove the 
frequency dependence, the equations for </ ? 5 nner (C) and ¥>2 nner (C) are simply the equations of 
motion for scalar fields in AdS2 with (mL 2 ) 2 = 0, 2 respectively. The radius of curvature of 
this effective AdS2 spacetime is related to the radius of curvature of the RN-AdS4 spacetime 
via L2 = L/y/6. 

It is simple to solve these inner region equations to give 



¥?r er (C) = a+ exp + ar exp f-^j , 

^(C)^4(c-0exp(^) + a 2 -(c + 0exp(-^ 



(46) 



where integration constants. The equations ( I45p are valid in the limit of small 

frequency and momentum (with respect to r Q /L 2 ) with (, defined in (jSJ), fixed. These 
requirements imply that r — tq < ro i.e. that these are near-horizon solutions. Ingoing 
boundary conditions at the horizon (£ — > 0) correspond to the choice a\ 2 = 0. At T = 0, 
Tq/L 2 is essentially the chemical potential fi of the field theory and so the following results 
are valid only for u,q fi. 



4.2 The outer region 

As in the T / case, we define the outer region by f[2"2"j) . Just as at T 7^ 0, the equations 
of motion in the outer region can then be solved to give the integral solutions (|30|) and (l3Tf , 
with Q = v3 for the T = case. In the limit T = 0, the fact that /(r) has a double pole 
means that these integrals can be done exactly but the results are lengthy and will not be 
presented here. To compute the Greens functions, we only require to know the behaviour of 
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the outer solutions in the matching region and near the boundary of the spacetime, and this 
will be presented in the following subsections. 



4.3 Matching 

To fix the integration constants c\ and c 2 in the outer region solutions, we need to know the 
inner and outer solutions in the matching region where the inner and outer regions overlap. 
The solutions overlap if we expand the inner solutions ( 146]) in the limit £ — > 00 ('far from 
the horizon') to obtain 



„mncr 
<Pl 



inner 
^2 



(4 + aT) [1 + (C 2 )] + ± (4 - af) [l + O (C 2 )] 



C(a+ + a 2 ) [i + o(r 2 )]+-^( 



648C 2 V ° 2 ° 2 

and the outer solutions fl30|) and fl3Tj) in the limit r — > r ('near the horizon' 



■) [1 + o (C 2 )] 



(47) 



outer 
Pi 



, outer 
¥>2 



(ci + 2c 2 ) u 2 
l8\/3r% (r — r ) 



(o) ,1 (o) . v / 3c 2 g 2 / 2 2 x 

Pi + ^2 + -^4- + { UJ C 1' W C 2) 



+ . . 



or (Ci - c 2 ) 
216rQ (r — ro) 



+ . . . + 



1 (°) , C2< i 2 , r>( 2 2 \ 



(r - r ) + 



(48) 

where we have only explicitly written the terms required for the matching. As was the case 
for u <C T, we are searching for a diffusion mode with u ~ q 2 and thus we can discard the 
O (lj 2 ci,(jJ 2 C2) terms as they are subleading at low q. 

After changing co-ordinates back from ( to r in (I4T1) . imposing ingoing boundary con- 
ditions fixes the ratio of the coefficients of the (r — r ) 1 and the (r — r )° terms of <^ 1 1 nner 
in the matching region, and the ratio of the coefficients of the (r — ro)~ 2 and the (r — ro) 1 
terms of y? 2 nner in the matching region as follows 



, „mncr „ H 



,„mner „H 
^2 = a 2 



1 + ^ (u 



+ ... 



r — r 
(r-r ) + ... + ^ R H 



(r - r )' 



(49) 



where we have absorbed a factor of uL 2 into and where 



~~ 6~~ 



3 r 6 



iu 6 L 
648 



(50) 



are (proportional to) the retarded Greens functions of scalar operators with conformal di- 
mensions A = 1, 2 in the CFTi dual to the near-horizon AdS2 geometry respectively. The 
relative coefficients of the terms in each series are determined by these CFTi Greens functions 
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simply because yj 1 1 nner and ip 2 nner behave like scalar fields in AdS2 with masses {mL 2 y = 0,2 
respectively. It is instructive for now to work with £{ R (u) and Q^(u), rather than using 
their explicit values ( 150]) . 

To determine C\ and c 2 for solutions which are ingoing at the horizon, we fix the coef- 
ficients of the power series in the outer solutions (jUJ) so that they agree with those of the 
inner solutions in the matching region (|49i) . Neglecting the O (w 2 Ci, w 2 c 2 ) terms, which are 
subleading at low q as previously explained, this fixes 



C-2 



Ayfal 2A^2,r u 2 ^ ) gf + £} R ( -MQf-tf<p { ? + r 3 W 2 y/3<p® + 3^ 



u 2 [120*V - r 2 (0{V + 2r u 2 ) 



12r 4 



l2Qf V f + rig™ ( V3^ } + <p 



,(°) 



12$ V - r 2 (0{V + 2r w 2 ) 



(51) 



Due to our previous approximations, these results are only valid to lowest order in u ~ q 2 



Recall that (u) is suppressed with respect to g[ (u) by a factor of u and therefore if 
we keep only the lowest order terms in these expressions (as we should for consistency), all 
dependence upon ^ R ( w ) drops out to give 



Ci = c 2 



6r 3 £{ R (u;) (vfy>£° ) + iff 



" 2 + 0? M 



(52) 



With these values for c\ and c 2 , the outer solutions correspond to fields which are ingoing 
at the horizon. 



4.4 Greens functions 



To determine the retarded Greens functions at T = 0, we substitute the outer solutions (1301) 
and (T5TT) (with Q = a/3 and C\ and c 2 given by (152"j) ) into the on-shell boundary action (TT9|) . 
This yields (where the argument (u) is shorthand for (u, q)) 



S 



j , - 



dudq 3rlq 2 g\ R 



>co (27r) 2 L 4 ^2 + |^£lR ( 



(-*) ^ ( W ) + ~fl> (- W ) M 



+ ^(ri 0) (--)^ 0) M + ^ 0) (--)^ 0) (-; 



(53) 
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up to contact terms. By using the relations (1411) and then applying the usual holographic 
prescriptions j56j[5T] , we find that the retarded Greens functions are 



R 3r V £ R (co) r 2 g 2 

Or 



JVJV ~*a q2L2y (54) 



12r 



r R _ rR _ (a;) _ v^rpg 2 

L^(^ + ^(o;)) 6* 2 (^-g 



and the Ward identities (117)1 are again trivially satisfied due to our gauge-invariant formu- 
lation. These results are valid to leading order at low to and q (with oj ~ g 2 ), and are valid 
only up to contact terms. 

The expressions fj54jl are precisely those obtained by taking the naive T = limit of 
the hydrodynamic results (T42]) ! In particular, they support a diffusion pole with dispersion 
relation given by (where the ellipses denote higher order terms in a small u, q expansion) 

u, 2 + Jh^ r (u;) + ... = — ► w = -,-^_ g 2 + ___ ; (55) 

which is simply the T = limit of the hydrodynamic relations ((!]) and 01]). As mentioned 
in the introduction, the main features of this dispersion relation are consistent with the 
principles apparent from previous studies of the low-energy physics of holographic theories 
|16 [[4"9"Ho"Tj . That is, the low-energy physics will be dominated by any gapless modes present 
in the theory in addition to the near-horizon geometry. The existence of such gapless modes 
is not apparent from the near-horizon geometry only but requires knowledge of the full 
bulk spacetime, although the decay rate of these modes is controlled by the near-horizon 
geometry. This is seen in the dispersion relation (|55p above - the existence of the diffusive 
mode required knowledge of the outer region of the geometry, but the power of momentum 
in the dispersion relation is controlled by the conformal dimension of the operator dual to 
the field (pi in the CFTx dual to the near-horizon AdS2 geometry. Unlike in the case of probe 
Dirac fermions studied in jl9], this dimension is fixed. 

It is instructive to compare our results to those of [37JI12], which analytically computed 
the leading behaviour of the retarded Greens functions in the limit u — > 0. Their calculation 
involves a matching procedure similar to that used here, but is formulated in terms of the 
gauge-invariant Kodama-Ishibashi fields - for which the equations of motion decouple [58] - 
rather than the fields <pi and (f2 used here. With these decoupled fields, one can determine the 
solutions in the near-horizon region more accurately than we have done here - in particular, 
the conformal dimensions of the CFTi operators dual to these fields in the near-horizon 
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region are given by 

1 1 a* I ~a% 

(56) 



1 1 

2 + 2 



\ 




These dimensions include q / /i-dependent corrections to the values 1 and 2 which we used in 
our calculation. The authors of [37] can therefore compute the uj — > scaling behaviour of 
the imaginary parts of the retarded Greens functions to find, for example, 

/ \ 2<S_-1 4 

T „ B I UJ\ uj ujq 1 , uj 

ImG^„ oc - = - + -*- log - + ..., 57 

where we have expanded in the limit q <C /i, compared to our result ImG!^^,, ocu/fi. 

Using the Kodama-Ishibashi variables also allows allows one to consistently include 
higher-order frequency corrections to the Greens functions which are analogous to the (uj) 
terms we neglected in our derivation. It is this contribution which gives a low-frequency dis- 
sipative conductivity (ignoring the delta function peak) of a {uj — > 0) ~ uj 2 [12] . Our results 
fl5^|) yield a (uj) = since Gj yJy (u, q = 0) vanishes due to our neglection of these higher-order 
terms. Despite these drawbacks of our approach as compared to using the Kodama-Ishibashi 
variables, the advantage is that we can solve the equations of motion in the outer region more 
precisely than has been previously done. As we have seen, it is this detailed knowledge of 
the fields far from the horizon which allows us to determine the existence of the diffusion 
mode. Of course, our results fl54l) give a 'zero temperature viscosity' equivalent to that found 
in [12] 

1 3r 2 1 s 

7] = - lim -lmG^ xyTxy (uj, q = 0) = — lim -Im^ R (uj) = —. (58) 

5 Greens functions at low, non-zero temperatures T < 

In this section, we will calculate the Greens functions (at leading order in small uj and q) at 
small (with respect to /i), non-zero temperatures which are outside of the usual hydrodynamic 
limit. In other words, in the range < T < uj fi. Heuristically, one would expect the 
calculation to proceed in a similar vein to the T = calculation but to now include small 
T/fj, and T/uj corrections. We will see that this is the case - for example the near-horizon 
AdS2 geometry is replaced by a near-horizon Schwarzschild-AdS2 geometry - and that the 
leading order Greens functions derived in this way are, in fact, equivalent to those obtained 
from the hydrodynamic calculation (H2|) . 
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5.1 The inner Schwarzschild-AdS2 region 



To define the inner region when T = 0, we wrote the equations in terms of the new radial 
co-ordinate C (defined in equation (jSj)) and then kept only the lowest order terms in uj 
and q in these equations. This corresponded to a near-horizon limit at small uj and q. To 
determine the leading T-dependent corrections to this result, we must keep T ~ uj in this 
limit also. To this end, we can replace Q (which determines T//i via equation (Q) with 



Co = 



r (3 - Q 2 ) AttT 



uj 



L 2 



uj 



(59) 



and keep this fixed in the low uj and q limit. Expanding linear combinations of the equations 
of motion f)18al) and f!18bl) at low uj and q, we find 



(0 + v\ (0 



12C + Co 



.C(6C + Co) 

(0 + ^; (0 [-..] + d (0 



+ 



+ ^ 2 (C) [...] + <P2 (C) [...]+¥>! (0 



1 



C 2 (6C + Co) 



+ . 



12C + Co 

c (ec + c ; 



+ . . . 



+ V2 (C) 



1 - 72C 2 - 12CCo 
C 2 (6C + Co) 2 



+ . . . 



(60) 

where the ellipses denote terms with positive powers of uj, q. As one might expect, in these 
limits the equations of motion decouple and are simply the equations of motion of scalars 
of masses 

(mL 2 ) = 0,2 in a Schwarzschild-AdS2 geometry respectively. This is most easily 
seen by changing co-ordinates to C = (p — Po) /ujL 2 with p = uL 2 ( /12, which yields the 
relevant Schwarzschild-AdS2 scalar equations of motion in a co-ordinate system in which the 
emblackening factor has the form f (p) = 1 — Po/p 2 - 
Solving these equations of motion, we find 



inner 



inner 
V?2 



(C) = at exp 
(0 = a 



-log (l + ^ 
Co H 6C 



> , Co i , 
C+ 12-6 ,eXP 



+ a x exp 



1 i M , Co 
T log 1 + — t 
Co V 6 C 



■flog(l + | 
Co V 6 C 



C + ^ + -jexp 
s 12 6 ' 



Co V 6 C 

(61) 



By an analogous argument to that at T = 0, these solutions are valid in the near-horizon 
region for small T ~ uj ~ q tq/ L 2 ~ Note that in the zero temperature limit Co — >■ 0, one 
can expand the logarithms and reproduce the T = results (1461) found previously. To impose 
ingoing boundary conditions, we can expand these inner solutions near the Schwarzschild- 
AdS2 horizon C — > to give 



, pinner 

(f l 2 ~ a 



1.2 



r 



'ro(3-« 3 ) 



n,2 



r 



-o(3-Q 2 ) 



(62) 



where the ~ indicates that we have absorbed various r-independent factors into the coeffi- 
cients af 2 - This is the expected result for fields at a horizon at non-zero temperature (recall 



21 



for example equation (|2"T|) ). Ingoing boundary conditions correspond to the choice a 12 = 0. 



We have already 



5.2 The outer region 

As in the previous cases, we define the outer region by the inequalities 
derived the solutions for ipl n 2 teT i n this region - they are given by equations (130 p and (13 ip . 

5.3 Matching 

Finding an overlapping region to match the inner and outer solutions is more involved in 
this case than in the previous cases and so we proceed carefully. Firstly, we expand the inner 
solutions f lBTj) near the boundary of the Schwarzschild-AdS2 spacetime (C — > oo) to give 



inner 
<Pl 



(C) = {at + al) + — (at- <h)[l + ... ] + ..., 



<4 nner (0 = C (4 + O + • • • + 



2592C 5 



[a+ (4 + 3< + Co 2 ) - <h ( 4 - 3<o + Co)] + • • • 



(63) 

where we have explicitly written only the leading terms whose ratio of coefficients is fixed 
by imposing ingoing boundary conditions at the horizon. After imposing ingoing boundary 
conditions at the horizon (a^ 2 = 0), absorbing constants into at 2 and changing back to the 
usual radial co-ordinate, we find that near the Schwarzschild-AdS2 boundary, 



, „ inner „H 

Vi = a i 



^ nner = <4 



IK 
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(64) 



where 



r>it 



IK 
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(65) 



are (proportional to) the retarded Greens functions of scalar operators with conformal dimen- 
sions A = 1,2 in the thermal CFTx dual to the near- horizon Schwarzschild-AdS2 geometry 
respectively. As in the T = case, the relative coefficients of the terms in each series are 
determined by these Greens functions simply because <£> 1 1 nner and (p 2 nner behave like scalar 
fields in Schwarzschild-AdS2 with masses {mL<i?' = 0, 2 respectively. 

Note that for a dimension 1 operator the Greens function is unchanged when one turns 
on the temperature T, whereas the Greens function of the dimension 2 operator receives non- 
zero T/oj corrections. These T/u corrections have the potential to introduce g-independent 
(but T-dependent) terms into the dispersion relation u (q) at low temperatures as we will 
shortly showjfl However, as in the T = case, the QJfip term is suppressed by a factor of to 2 



J For a toy model of probe D-brane systems similar to those of [251153], it is similar T/lj corrections to the 
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with respect to the Q^t contribution (recall that we assume u > T in this section) and thus 
makes only a subleading contribution to the dispersion relation. It is for this reason that 
the dispersion relation does not acquire any significant g-independent (but T-dependent) 
contributions. 

To show this, we must match these inner solutions to the outer solutions over a range of 
r where both solutions are valid. Simply expanding the outer solutions around r = tq yields 
the results (1321) and (1331) which do not overlap with the inner solutions (1641) . This is because, 



as discussed in section I3.3[ these solutions are valid in the usual hydrodynamic limit w<gT 
whereas the inner solutions (164"|) are valid for u > T. The correct procedure is to expand 
the outer solutions in the limits r - r « r and 3 — Q 2 <C 1, where this second limit is 
equivalent to T//i < 1, which is also the range of validity of the inner solutions ( 16~4"|) . In 
these limits, 

f(r) — ► 1 (r - r ) 2 (r 2 + 2rr + 3r 2 ) = f (r) , (66) 
where the subscript indicates that f(r) takes on its T = form, and 

_ 4Q 2 r 4Q 2 / _ rp\ 

3(1 + Q 2 )r 3(1 + Q 2 )V r)' 1 ' 



Note that in these limits, the outer region inequality (122]) no longer implies that uj <C T. 

With these substitutions, the outer solutions f[3"Uj) and f[3"Tj) can be expanded near the 
horizon r = tq, where they overlap with the inner matching solutions (IMl) just derived. In 
this matching region, the outer solutions take the form 



, outer \ - 

^2 = - — 72 + • • • + 

216rQ (r — ro) 
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uj 2 (ci - c 2 ) 
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(68) 

where, as before, we have only written explicitly the terms which are required for matching, 
and have neglected the O (w 2 Ci, cj 2 C2) terms as they are subleading in a low q expansion with 
u ~ q 2 . 

To impose incoming boundary conditions at the horizon on the outer solutions (130 j) and 
( 13 ip . we simply compare the expansions above to those of the ingoing inner solutions in the 
matching region (IM]) . This fixes c\ and C2 to be (at lowest order in an expansion in u and q) 



6Qr 3 G\ R T (u 
ci = c 2 = 



(0) 4Q (0) 

^1 + 3(1+Q 2 )^2 



(69) 



inner region solution which result in the Fermi liquid-like zero sound dispersion relation — Im (cj) ^ k 2 + T 2 
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As in the T = case, C/^r ( u ) contributes to these constants only at a subleading order in 
the low uj and q expansion, and is thus outside the range of validity of our approach. 



5.4 Greens functions 

Our final step in computing the retarded Greens functions is to substitute the outer solutions 
f )30|) and f l3Tj) with ingoing boundary conditions given by fl69|) into the on-shell action ffl~9jh 
This gives (where the argument (u) is shorthand for (u, q)) 



. , dudq 3r 2 q 2 G\ R T (uj 



(o) / n (o) / n . 4Q / (o) , \ (o) / \ i (o) / \ (o) / A 
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(70) 

up to contact terms. After using equation (j4Tj) to express the boundary values of the gauge- 
invariant fields in terms of the sources of dual field theory operators, we can compute the 
retarded Greens functions of these dual operators via the usual prescription [56|l57]. The 
results, valid for T < u <C /i, are identical to those in the hydrodynamic limit given in 
equation (l4"2"j) and the Ward identities (IT7|) are automatically satisfied due to our gauge- 
invariant formulation. These results are valid only at lowest order in T, u and q (with 
T, uj ~ q 2 ) and are valid only up to contact terms. 

In particular we note that when T < u <C //, there is a diffusive mode with dispersion 
relation (where the ellipses denote terms which are higher order in q and T) 

w2 + Mr?OT 6 " M + - = - w = - ! 3mtW 2 + -- (71) 

which agrees precisely with the diffusion mode predicted by the hydrodynamic derivative 
expansion in equations (GQ) and (jl]). As in the T = case, it is the dimension of the CFTi 
operator dual to the field tpi in the near-horizon Schwarzschild-AdS2 geometry which controls 
the power of q in the dispersion relation. As we indicated earlier, it is the lack of any In- 
dependence of the CFTi Greens function Q]^ (u) which result in there being no terms in 
the leading order dispersion relation (ITTj) which are q- independent but T-dependent. 



6 Numerical computations of poles and spectral func- 
tions 

As a check of the validity of our analytic results, we have computed numerically both the 
poles of the retarded Greens functions and the spectral functions of the operators studied 
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above. This also allows us to reconcile our results with some previous numerical results for 
the planar RN-ACIS4 theory [37J[39]. We have used two different numerical methods to do 
these calculations. 

The simplest method conceptually is to numerically solve the equations of motion ( 1 18 aft 
and f!18bj) by integrating out from the horizon of the planar RN-AdS4 spacetime. If one 
demands ingoing boundary conditions at the starting point of the integration (the horizon), 
then both the poles of the retarded Greens functions and the spectral functions are easy to 
determine from the values of the fields at the end point of the integration (the boundary). 
This numerical procedure is described in detail in |57j . 

It is well-known that this method becomes very difficult to implement when T = 0, as 
the horizon becomes an irregular singular point of the equations of motion (when T 7^ it 
is a regular singular point). Leaver's method [60] is an alternative way to compute the poles 
of the Green's functions which does not require integration around a singular point. This 
method involves making a quasinormal mode ansatz for the fields <p (r) = <^qnm ( r ) Wm. (r), 
where v?qnm (r) contains the relevant asymptotic boundary conditions such that if (r) is 
ingoing at the horizon and its leading term vanishes at the boundary of the spacetime. After 
making this ansatz, one then expands the remainder ip rem , (r) as a power series (up to order 
-/Vmax.) around a regular point of the differential equations - usually taken to be the midpoint 
between the horizon and the boundary f| One can then expand the equations of motion in a 
power series up to order iV max around the same point to give iV max . + 1 algebraic equations for 
the iV max + 1 coefficients of the power series expansion of <f rcm , (r). If a non-trivial solution 
to this set of equations exists, then it implies the existence of a bulk quasinormal mode and 
hence a pole in the retarded Green's functions of the dual field theory. Note that although 
this method is most useful at T = where the method of direct numerical integration is 
difficult, it can also be applied at non-zero T. This numerical procedure is described in more 
detail in 



6.1 Zero temperature poles 

As indicated above, we will use Leaver's method to determine the poles of the Greens func- 
tions when T = 0. For the planar RN-AdS4 theory, it is simplest to implement this method 
using the gauge-invariant Kodama-Ishibashi variables (rather than the variables ffTBT) ) for 
which the equations of motion decouple [58] • This was done in detail in [37] and we refer the 
reader to that work for the technical details, as our calculation is identical to that described 
therein. 

Leaver's method indicates the existence of a discrete spectrum of poles with complex fre- 
quencies, as well as a very large number of closely-spaced poles along the negative imaginary 

4 For this method to work, the closest singularities to this regular point should be those at the horizon 
and the boundary of the spacetime. 
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Figure 1: The locations of the purely imaginary 'poles' nearest the origin as determined by 
Leaver's method when q/n = 0.1 (shown as dots), as a function of the order of the power 
series expansion used iV max .. The red line shows the analytic result ( 153|) for the T = 
diffusion mode. 

frequency axis. As described in [37], these purely imaginary 'poles' are actually a numerical 
representation of the fact that there is a branch cut along the negative imaginary u axis due 
to multi- valued terms in the Greens functions (which can be seen in equation (157|) ). If one 
increases the accuracy of the implementation of Leaver's method by increasing iV max ., the 
spectrum of these imaginary 'poles' becomes denser and denser as they move closer to the 
origin and in the limit iV max . — > oo, this spectrum of 'poles' will become a continuous line 
along the negative u axis, which represents the branch cut. 

In figure [U we show how the locations of the purely imaginary 'poles' are dependent 
upon the accuracy of the calculation iV max .. As expected, most of these 'poles' have a 
location which is not stable - as one increases iV max ., they move closer together. However, 
the location of the pole closest to the origin is stable as one increases iV max . and its location 
agrees precisely with that of the zero temperature diffusion mode predicted by our analytic 
calculation f )55|) . We thus conclude that in addition to the branch cut, there is a genuine 
diffusive pole of the retarded Greens functions with dispersion relation flB3|) . Our analytic 
calculation has captured this diffusive pole but not the multi-valued term (which is higher 
order in u,q). In this sense, the analytic calculation of the low frequency Greens functions 
in [37] is complementary to ours, as it determines the existence of the multi-valued term, 
but not of the diffusive pole. 

In figure El we show the dispersion relations of some of the purely imaginary 'poles' at 
fixed iV max . = 300. The location of the 'poles' which are just a representation of the branch 
cut are approximately g-independent. However, the genuine pole has a dispersion relation 
which is described accurately by our analytic diffusion result f )55|) when q -C [A. 
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Figure 2: Dispersion relations of the purely imaginary 'poles' nearest the origin as determined 
by Leaver's method, at fixed iV maXi = 300 (shown as dots). The red line shows the analytic 
result fl55|) for the T = diffusion mode. The data on both plots is the same - the left hand 
plot is a magnified version of the near-origin region of the right hand plot. 



6.2 Non-zero temperature poles 

At non-zero temperature, where the branch cut is resolved into a discrete spectrum of poles, 
it is possible to compute the spectrum of poles using either Leaver's method (as was done in 
[37|f39]) or direct numerical integration. We computed numerically the dispersion relation of 
the candidate diffusion mode and fitted this result to a function of the form u = —iVq 2 —iAq 4 , 
for various values of T/fi. The corresponding numerical values of the diffusion constant T> as 
a function of T j /i are shown in figure [3j The results obtained using each numerical method 
coincide, and are well-described for all temperatures by the hydrodynamic results ([T]) and 
(jl]). These results differ from those in figure 9 of [39] - this is because the fitting function 
used there does not include a q 4 term, the effects of which are important towards the upper 
range of the values of q/fi that we have studied numerically. 



6.3 Spectral functions 

As a further check of our analytic calculations, we have computed numerically the T ^ 
spectral functions of the transverse components of and J M , which are defined as 

Xoo (w, q) = -2lmGao (w, q) ■ (72) 



This was done by directly integrating the equations of motion f)18ap and f!18bj) from the hori- 



zon to the boundary as described in [57] . The resulting spectral functions, at all numerically- 
accessible temperatures, are well-described by the spectral function formulae computed from 
our analytic results (H2l) . In figure H] we show a sample of these results although for concise- 
ness we will not give a detailed exhibition of them here. 
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Figure 3: The numerical values of the dimensionless diffusion constant [iD (shown as dots) 
as a function of T//x as computed by iV max . = 300 Leaver's method (left hand plot) and 
direct numerical integration from the horizon (right hand plot). These were extracted by 
fitting the function u> (q) = —iDq 2 — iAq 4 to the numerical location of the pole u (q) at 
qj 'fi = {0.01, 0.05, 0.1, 0.2, 0.3, 0.4, 0.5}. The red lines show the hydrodynamic result for jiD, 
given by equations ([T]) and (HI, which according to our previous analysis should be valid for 
all T//i shown. 




Figure 4: Numerical results for the = 0.1 spectral functions xt*«t*» (shown as dots) and 
Xjyjy (shown as crosses), computed via direct numerical integration from the horizon, for 
T I jj, m 0.005 (left hand plot) and T//i fa 0.16 (right hand plot). Xt*vt*v is plotted in units 
of 7"o/2k4L 4 and Xjvjv is plotted in units of Tq/2k\. The analytic results (T4"2"|) for xr*»r*» 
(dashed line) and Xjvjv (solid line) are plotted in red. 
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7 Discussion 



In this paper we have studied the two-point functions of the transverse components of T^ u 
and J M for the field theory dual to the planar RN-AdS4 black hole. Our main result is 
that a long-lived diffusion mode - described by the hydrodynamic formulae (pQ) and (jlj) - 
exists outside of the commonly-studied range u <^ T. Specifically, it exists provided that 
uj <T or \i. It exists even in the extremal limit T = 0. 

Recall that hydrodynamics is constructed as a derivative expansion, valid when the length 
scale over which macroscopic quantities vary is much larger than a characteristic length scale 
/ m f p . Our results indicate that this length scale is given by 

I = S = ^ = 1 Q (73) 
mfp 4vr( e + P) 3r (l + Q2) 3^(1 + Q 2 )' 

since with this definition the hydrodynamic diffusion result is valid when cuZ m f p , ql m f p <C 1. 
Note that Q is implicitly a complicated function of T//x, given in equation flU]). In the \i = 
limit, Z m f p ~ 1/T and in the T = limit, Z m f p ~ (we have neglected numerical factors of 
order 1). By definition, this mean free path is equivalent to that suggested in [20] for related 
global black hole solutions. Our work confirms the assumption therein that the lifetime of 
the longest-lived mode is proportional to the quantity s/(e + P) at all temperatures. 

A natural extension of the work we have done here is to calculate the low-energy ex- 
citations of the longitudinal components of T^ v and J M for this theory outside the usual 
hydrodynamic range u T. We expect that the results will be similar in spirit to those 
presented here in that there will be a sound mode with dispersion relation given by (j2J) and 
(jlj), and a charge diffusion mode, provided that u) <C T or \l. There is already numerical 
evidence to support this claim - in [38J it was shown numerically that there is a T = sound 
mode whose real part agrees closely with the hydrodynamic result ([2]) and whose imaginary 
part differs from the hydrodynamic result by around 10%@ Furthermore, in [10] the sound 
dispersion mode was studied at non-zero temperatures. Although not explicitly stated there, 
both the real and imaginary parts of this mode are consistent with the hydrodynamic results 
([2]) and (jlj) for all T provided that A convincing verification of this expectation will 

require an analytic derivation of the dispersion relation of the low-energy modes, and we are 
hopeful that the method we have used here may be suitably modified to provide this. 

The low-temperature two-point functions that we have derived are valid only at lowest 
order in an expansion in small u,q and T with u,T ~ q 2 . It would be interesting to try 
to generalise our low-temperature calculations to include higher order terms in u, q and T. 



5 No T = charge diffusion mode was identified, but this may just be due to the difficulties in numerically 
disentangling it from the branch cut as described in section [SJ 

6 This can be seen by comparing the numerical results for the sound attenuation in hgure 1 1 of [3D] to 
those for the diffusion constant in hgure[3]of this paper, and noting the relation between these two quantitites 
predicted by hydrodynamics from equations (fTJ) and ([2]). 
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Firstly, this will allow us to determine more accurately many of the 00 — > properties of our 
theory such as the conductivity and the scaling of the imaginary parts of the Greens functions 
in this limit. Presently, the calculations of [371112] are more accurate in this regard. Perhaps 
more importantly, this would indicate whether the predictions of hydrodynamics are valid 
at low temperatures to higher orders in the derivative expansion, or just at the leading order 
as shown here. Due to the presence of the branch cut, we expect that hydrodynamics will 
in fact break down at some higher order. In a similar vein, one could investigate the range 
of validity of non-linear hydrodynamics in this theory by computing higher-point Greens 
functions. 

Another obvious question is whether the applicability of hydrodynamics when T < u <C 
fi, as shown in this paper, is generic to (non-probe D-brane) holographic field theory states 
with a large chemical potential. The planar RN-AdS4 solution is in some ways a very special 
solution, as its entropy density s is non-zero when T = 0. This feature is intimately linked 
to its near-horizon AdS2 x M. 2 geometry which is critical for our results since the decay rate 
of the diffusion mode is controlled by the dimension of a scalar operator in the CFTi dual to 
the AdS2 geometry. States with s = at T = can be studied by including matter fields in 
the gravitational theory - these states typically have a different near-horizon geometry which 
can, for example, have Lifshitz scaling and/or violate hyperscaling. A particularly interesting 
example are the near-horizon geometries of (62] which are conformal to AdS2 x M. D ~ l and 
have a vanishing zero temperature entropy. 

It would be of interest to determine which (if any) of the hydrodynamic results ([1]) and (j3J) 
break down at small T for different near-horizon geometries, and if the applicability of these 
formulae at small T is dependent upon the extent of fractionalisation of the U(l) charge in 
the dual field theory (which has a strong influence on the near-horizon geometry as reviewed 
in [63]). It would be a surprising result if the hydrodynamic results (JTJ, (T5]) and (j2J) are valid 
at low temperatures irrespective of the dual geometry. This would mean that the leading q 2 
term in the decay rates of the diffusion and sound modes in 'realistic' field theories (i.e. those 
with s = at T = 0) vanish at T = and hence these modes would live for an anomalously 
long time compared to, for example, the well-understood zero temperature sound mode of 
a Fermi liquid [2"T j l2"3" | l2"l] . However, as we have stressed previously, one of the main lessons 
learned so far from studying strongly-coupled field theories using gauge/gravity duality is 
that their transport properties can be very different from those one would expect from a 
quasiparticle-based theory. 
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